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ABSTRACT 

We use data from the Sloan Digital Sky Survey and visual classifications of morphology from the 
Galaxy Zoo project to study black hole growth in the nearby Universe {z < 0.05) and to break 
down the AGN host galaxy population by color, stellar mass and morphology. We find that black 
hole growth at luminosities L[Oiii] > 10 ergs~^ in early- and late- type galaxies is fundamentally 
different. AGN host galaxies as a population have a broad range of stellar masses (lO^*' — IO^^Mq), 
reside in the green valley of the color-mass diagram and their central black holes have median masses 
around lO^'^M©. However, by comparing early- and late-type AGN host galaxies to their non-active 
counterparts, we find several key differences: in early-type galaxies, it is preferentially the galaxies 
with the least massive black holes that are growing, while late-type galaxies, it is preferentially the 
most massive black holes that are growing. The duty cycle of AGN in early-type galaxies is strongly 
peaked in the green valley below the low-mass end (lO^^M©) of the red sequence at stellar masses 
where there is a steady supply of blue cloud progenitors. The duty cycle of AGN in late- type galaxies 
on the other hand peaks in massive (lO^^M©) green and red late-types which generally do not have 
a corresponding blue cloud population of similar mass. At high Eddington ratios {L/L-^dd > 0.1), 
the only population with a substantial fraction of AGN are the low-mass green valley early-type 
galaxies. Finally, the Milky Way likely resides in the "sweet spot" on the color-mass diagram where 
the AGN duty cycle of late- type galaxies is highest. We discuss the implications of these results for 
our understanding of the role of AGN in the evolution of galaxies. 

Subject headings: galaxies: evolution; galaxies: formation galaxies: Seyfert; galaxies: active 
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1. INTRODUCTION 



The nature of the co-evolution between galaxies and 
their supermassive black holes is still a mystery a 
decade after the discovery of the rela tionship between 
spheroid mass and b l ack hole mass ([Magorrian et alj 
l998t iGebhardt et all [2OOOI: iFer rarese fc Merrittl 120001: 



Tremaine et alJ 120021 : iHaring fc Rix 2003). Even more 
mysterious is its connection to the phenomenon of 
downsizing, the anti-hierarchical growth of galaxies and 
black holes. While dark matter grows and assembles in 
a strictly hierarchical fashion, with the most massive 
dark matter haloes continuing to grow over cosmic 
time by accretion of smaller haloes, the formation and 
growth of galaxies inside those halos proceeds differently 
e.g rSilki 119771: iRees fc Ostrikcr 1977; White fc Reel 
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Fall fc Efstathiou .1980; Blumenthal et al., 1984 
Baryonic matter is converted 



Cole fc Lac ev 1996 



into stars within dark matter haoles and regulated by 
feedback processes. The stars in the most massive 
galaxies were not only formed first, but also on the 
shortest time scales, with further star formation at 
later epochs being suppressed. Less massive galax- 
ies have increasingly more extended star formation 
histories as star formation in them, and their progen- 
itors, proceeded over increasingly longer timescales. 
This picture has been verified both by "archeological 
studies" of local galaxies, and by prob ing ongoing 
star fo rmation at in cr easing redshifts (ICowie et al.l 
1991 [Heavens et al.l [200I iThomas et all 120051: 



20061 Jimenez et al 



Nelan et al.l 120051: iBern ardi et al.l 1200 61: iBundy et al l 



y et ai.l 

.2007: Perez-Gon zalez et alTT2008l: 
Collins et al.l 120091: iTho'mas et al.ll2009t ). and has alter- 
natively been called 'downsizing' or 'anti-hierarchical' 
growth. 

The growth of supermassive black holes seems to 
broadly mirror the anti-hierarchical trend of their host 
galaxies. Quasars powered by ~ IO^Mq black holes are 
already observed about 1 Gyr after the Big Bang, im- 
plying that they must have reached t hese masses in a 
very brief, intense gr owth phase (e.g., iFan et all I2001L 
[20031 [2OOI [Fan[[20Ci6[ ). Studies of the luminosity func- 
tion of optical- and X-ray-selected active galactic nu- 
clei (AGN) in deep fields observed by Chandra and 
XMM-Newton have revealed that the most luminous 
AGN were most common in the early Universe, while 



the lower-lunrinosity AGN peak l ater (lUeda et al 



Merlonill200l iBarger et al.H2005l: iLa Franca eta! 
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Croom et al.|[2009[ ). In other words, the pattern of black 



hole growth inferred from the luminosity evolution of 
AGN seems to exhibit a qualitatively similar trend as 
galaxy growth. Studies attempting to rec over the history 
of black hole growth supp ort this picture ([Marconi et al.l 
[200I IShankar et al.l[2009h . 

It has been suggested in recent years that energy in- 
put from AGN not only self-regulates the growth of both 
black holes and their host galaxies (I S ilk fc ReesI 119981 : 
[Kind I2003D to give rise to the observed Mbh — cr re - 
lation ([Gebhardt et al.ir2000l: IFerrarese fc Merrittl [2OOOI) 
but that this AGN feedback also, as [Bower et al. 
described it, "breaks the hierarchy of galaxy formation" 
giving rise to the observed anti-hierarchical nature of 
galaxy star formation histories. However, the precise na- 
ture of this connection is not understood, and it is not 



clear if all or only some AGN phases affect their host 
galaxies. 

In this paper, we explore properties of the host galax- 
ies of AGN and the nature of their black hole growth 
in the local Universe {z < 0.05) and discuss what these 
properties imply for the role of AGN in the evolution 
of the host galaxies. In our investigation, we go beyond 
characterizing AGN host gala xies as a single class, as 
in some previous stu dies (e.g.. iKauffmann et al.|[2003a[ : 
iHeckman et al.l [20041 ). and dissect the AGN host galaxy 
population into sub-populations by observed properties 
such as mass, color and most importantl y, morphology 
(e.g., [Best et al.|[2005[: iHickox et al][2009[ ). We use the 
visual classifications of morphology by cit izen scientists 
taking part in the Galaxy Zoo project ([Lintott et al.l 
[20081 to divide both AGN host galaxies and their quies- 
cent counterparts into morphology classes to reveal sev- 
eral new and surprising properties of AGN host galaxies 
that have implications for our picture of the galaxy-black 
hole co-evolution. Our main conclusion is that the early- 
and late-type host galaxies of AGN are, despite appar- 
ent similarities, very different populations, and that the 
AGN phases occurring in them play a very different role 
in their co-evolution. 

In Section [21 we describe our sample selection and the 
derivation of galaxy properties such as stellar masses and 
morphologies. In Section[3l we motivate and describe the 
selection of AGN via emission line diagnostic diagrams 
and the derivation of accretion properties and black hole 
masses. We present our observational results in Section 
[31 discuss their implications in Section [5l and summarize 
in Section [HI 

Throughout this work, we assume cosmological param- 
eters (rim — 0.3, r^A = 0.7, Hq — 70), consistent with 
the Wilkinson Microwave Anisotropy Pro be Third Year 
results in combination with other data ([Spergel et al.l 
[2OOI . 

2. DATA AND SAMPLE SELECTION 

In this section, we describe the selection of both the 
parent galaxy population and the AGN host galaxy sam- 
ple we use in this paper. 

2.1. SDSS Photometric and Spectroscopic Data 

The photometric and spectroscopic data used here are 
taken from the S l oan Digital Sky Survey (SDSS) DR7 
(lYork et al.l 120001 : iStrauss et all l2002t [Abazaiian.eLalJ 
I20090 . Throughout the paper, we use modelMags for 
galaxy colors, as they are derived from the best-fitting 
exponential or de Vaucouleurs galaxy profile and are thus 
best suited for measuring galaxy colors. 

2.2. Galaxy Sample Selection Criteria 

Our sample selection criteria ensure a large and com- 
plete sample of galaxies to low stellar masses whose SDSS 
images are suitable for detailed morphological classifica- 
tion. We select all galaxies with SDSS spectra classified 
as GALAXY (Stra uss et al. 2002 ) in the redshift interval 
0.02 < z < 0.05 and limit our selection to galaxies with 
r < 17 AB mag. In order to create a magnitude-limited 
sample that is close to a stellar mass-limited sample, we 
further limit this sample to Mz^petro < —19.5 AB mag. 
We choose the z-band rather than the r-band since the 



Galaxy Zoo: AGN Host Galaxies 



3 





5877^t^4(l1S9^MJ73 




■ * 

■.".••vW..- 


7 

5877i6«1(»I71l7«3 


5877356H81!OJ!i1B 


D 


* 

% 








»S7TaH>1sa7t77Sj16 



58771*983902 








I. 

5877I5(l7«9l»»aM 


* 

5877i6«05i59BJ+90 


I. m ^ 

Mn73967ew4«aoe4ir 
















3877»1«Hl30iJ*7t 


SS7rtWJ44l*SsSS01 




ai7Ti71»6 163681 33 




1 ^ 

SS7rt«l7i7«73l4lS 



5877Ii9844S((«9ie I 5877^t^4(l15S«^W35 I 5a373W79+1J10917« I 5877I7Ji1W«OS»79 I 56fra794S9e0«1t-*0a I WUTiaiHWBaASOW: 





3S7i'i)*i1l>»*»01rrj3 . 



Fig. 1. — Example images of AGN host galaxies in three morphology classes. The top two rows arc late-type host galaxies, the middle 
rows are indeterminate-type and the bottom two rows are early-typc host galaxies. In the top left of each image, we show a bar of 5 kpc 
length at the redshift of the object. In the bottom, we provide the SDSS object ID. 



redder z-band is closer to stellar mass. This yields a total 
of 47,675 galaxies. 

T he SDSS fibres are circular with a diameter of 
3" (|York et al.|[2000t IStrauss et al.|[200l . which corre- 
sponds to a physical footprint of 1.2 - 2.9 kpc in diame- 
ter at the lower and upper redshift limits of our sample. 
The spectroscopic information used in this work therefore 
samples only the central regions of galaxies, rather than 



the bulk light which is sampled in more distant (z 0.1) 
galaxies. We describe the selection of AGN from this 
parent sample in Section |3l 

2.3. Galaxy Zoo Morphologies 

In the past, obtaining accurate morphologies for large 
samples of galaxies has been difficult. Only with the in- 
volvement of over 200,000 citizen scientists in the Galaxy 
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TABLE 1 

Galaxy morphology statistics from Galaxy 
Zoo 



Galaxy morpholgy class 


Number 


Percentage 


Late-type 


16246 


34.1 


Indeterminate-type 


22483 


47.2 


Early-type 


8928 


18.7 


All 


47675 


100 



Zoo project (' Lintott et alj [20081 ) have we been able to 
gather muhiple independent visual classifications for the 
entire SDSS main galaxy sample that are comparable to 
those of professional astronomers. Galaxy Zoo provides 
a web interfac4H3 where volunteers from the general pub- 
lic are shown randomly selected images from the SDSS 
main galaxy sample. Since every object supplied for in- 
spection by Galaxy Zoo has been viewed multiple times 
by independent classifiers, we are able to use the distri- 
bution of votes to choose the purity of the classifications 
of individual galaxies by specifying the required level of 
agreement among classifiers. 

We use t he d efinition for the clean sample of 
iLand et al.l ()2008[ ). which requires an 80% majority 
agreement on the morphology of any object. We com- 
bine all spirals regardless of orientation into the late- 
type class and all the spheroids into the early-type class, 
while those objects without a consensus at the 80% level 
are placed in the indeterminate class. We use the label 
'early- type' rather than 'elliptical', as we specifically in- 
clude lenticular galaxies in this class, while the late-type 
class contains anything with a discernible disk, from an 
Sa to an Sd galaxy. 

We include merging galaxies with the indeterminate- 
type class; in any cas e, major mergers in the local Uni- 
verse are rare (1-2%; iDarg et al.ll2009U and references 
therein) , and the incidence of emission- line AGN in them 
is not significantly enhanced over the general population 
(|Darg et al.ll2009a[ ). In Table [H we present the total 
numbers and population fractions of the various mor- 
phology classes defined by the Galaxy Zoo classifications 

In order to give a sense what these three morphological 
classes mean, we show in Figure [T] a randomly selected 
sample of AGN host galaxies in each morphology class 
class. We stress that these are AGN host galaxies and 
not objects taken from the entire population. For exam- 
ples of image s of no rmal galaxies, we refer the reader to 
iLintott et atl H2008D . 

2.4. Stellar Mass Measurements 

We measure stellar masses for all objects in this sam- 
ple by fitting the five SDSS photometric bands to a li- 
brary of 6.8 X 10^ model star form ation histories gen- 
erated from iMarastonI (|1998l , I2005D stellar models. The 
star formation histories are parameterized as tw o bursts 
(e.g.. lSchawinski et all2007tlKavirai et al.ll2007[ ) of vary- 
ing ages and mass-fractions, fixed solar m etallicity, vary- 
ing d ust extinction (0 <E(i? - V)< 0.5; iCalzetti et al.l 
120001 ) and exponential e-fol ding time (r = 10, 100, 1000 
Myr); all models assume a iSal peteil (|1955[) initial mass 

■^^ See http : // zool.galaxyzoo.org 



function. Stellar masses are measured by finding the min- 
imum of the statistic in the parameter space probed, 
while the errors are determined from the range of stellar 
masses within 1 — cr of the minimum. We do not vary 
the initial mass function. 

We test the reliability of our stellar mass measurements 
for AGN host galaxies with the methodology described 
in Appendix [X] We find that our stellar mass measure- 
ments are robust for the AGN present in our sample. 

3. AGN SELECTION 

The selection and characterization of AGN is a chal- 
lenging problem, as no single wavelength regime al- 
lows for a perfectly unbiased sample. Various mid- 
infrared selection tec hniques have been proposed based 
on Spitzer data ('e.g..'L acv et al.l 20041: Stern et al.ll2005l : 
Donley ct al. 2005; Hic kox et al. |200£ ) . but these can 
be inefficient depending on depth, missing a large frac- 
tion of genuine A GN while selecting a substantial num- 
ber of non-AGN (jCardamone et all [2001 ) . Hard X-ray 
selection {E > 10 keV) is the least biased selection 
technique, as it can detect even the most heavily ob- 
scured, Compton-thick systems whose existence is re- 
quired by the X-ray background (jTreister fc Urrvll2005l : 
iTreister et al.l[2009al ). However, current surveys of the 
hard X-ray sky remain too shallow or are confined to 
small areas compar ed to SDSS (e g., iTueller et al.l [200l 
[Winter et al.ll2008l: [^Reister et al.|[2009air 

For this reason, we resort to a less complete, but 
more sensitive selection technique than hard X-rays 
by searching for galaxies exhibiting narrow emission 
lines originating from a cone of radiation photoionize d 
by the centr al engine (e.g.. lEvans et al.l I1991L Il993| ). 



First used by iBaldwin. Phillips, fc TerlevichI (|198lD and 
extended by 'Veil leux fc OsterbrocU ( 1987D . emission 
line diagrams have been improved and refined in their 
interpretation b y many theorists using photoionisa- 
tion codes (e.g., iStasinskal Il990t iLeitherer fc Heckma" 



19951 ILeithere^r et al l 119991 : IKewlev et al.l I2001L 1200 



Groves et al.l I2 004a"t7) and by observational work 



(e.g.. [Veilleux fc Ostcrbrock 1987; Ho et al.l 



i997a''bl'c': 'Kauffma nn et al.ll2 003a: St asihska et al.l 



1995 



2006 



r 

iSarzi et al..,2006; Stasihska et al... 2008( 1 

The emission line selection we use in this work does 
not include unobscured (Type 1) AGN with broad lines 
and may be biased against the most highly obscured 
AGN. If the Unified M odel of AGN ([Antonuccil [1991 
lUrrv fc Padovanill995( ) is correct and these difference are 
due to the orientation of the putative obscuring medium 
in the central engine, and if the obscuration does not 
cover the narrow-line region, then narrow-line (Type 2) 
AGN are a fair sub-set of the entire AGN population 
within the luminosity range covered. 

We now discuss in detail how we measure emission 
lines, use them on emission line diagrams, and assess 
the completeness of this selection. The final result of our 
selection described in this Section is that of the high- 
luminosity end of the local AGN population, and there- 
fore of sites of substantial black hole growth. For rea- 
sons of completeness, a similar analysis of low-luminosity 
AGN with significantly sub-Eddington accretion rates is 
not possible via emission line selection due to contami- 
nation from star formation (see Section [3731 for details). 
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3.1. Emission Line Measurement 

We measure the emission line fluxes using an anal- 
ysis tool called Gas A ND Absorpt i on Li ne Fitting al- 
gorithm (GANDALF; iSarzi et all I2006D . which uses 
the Penalized Pixel-Fit t ing m ethod PPXF method of 
iCappellari fc EmsellemI ()2004D to simultaneously fit the 
stellar continuum and nebular emission 

linetEE- The ad- 
vantage of this tool is that it is designed to minimize 
template mismatch, where e mission and absorption lines 
are difficult to separ ate (see iFalcon-Barroso et al.l 120061 
and ISarzi et al.l[2006l for exampl e applications). We use 
a set of stellar templates from iTremonti et al.l ()2004l ) 
plus Gaussian emission line templates. We measure 
[Oiii] A5007 fluxes for each processed SDSS spectrum 
and compute the total luminosity L[Oiii] by correcting 
for extinction based on the measured Balmer decrement. 
GANDALF also provides measurements of the stellar 
velocity dispersion. Since the fixed SDSS spectroscopic 
fibre samples the light out to a different fraction of the 
eff ective radius, we use th e empirical correction derived 
bv ICappellari et aT] ()2006D to correct the velocity disper- 
sion to the effective radius of the bulge surface-brightness 
profile. 

3.2. Emission Line Diagnostics and AGN Selection 

Our AGN selectio n technique is the same as used 
by iSchawinski et al.l ((2007) : we first separate galaxies 
whose narrow emission lines are dominated partially 
or wholly by star formation on the [Nil] /Ha diagram 
(Figure |2l lef t ) by using the extreme starburst line of 
iKewley et al.l (|2001[ ) to retain only those beyond this 
line (i.e., the top-right of the diagram). Between this 
extreme starb urst line an d the empirical pure starburst 
line of Kauffm ann et al.| ((lOOSb), there is a substantial 
population of composite objects where both AGN and 
star formation are comparable in ionizing luminosity. We 
exclude this class for this study and note that by remov- 
ing this class of potential AGN, we may be removing an 
important phase in the AGN-galaxy co-evolution at low 
redshift. The main problem with the composite class is 
the lack of easily accessible spectral indicators that would 
allow a clean separation of AGN and star formation. 

We are thus left with galaxies whose emission lines 
are dominated by sources of ionization other than young 
stars. These sources may be AGN, but can also include 
slow-moving shocks and gas excited by post-AGB stars 
and horizontal branch s tars (sec Ho 2008; Stasiiiska et al.l 
I200I ISarzi eraI1l2009l) . These non-stellar sources are 
empirically divided into two branches, a division that is 
most obvious in the [OiJ/Ha diagram (Figure [51 right). 
The lower branch of this population was first designated 
as 'low-ionization narrow emission-line regions', or LIN- 
ERs, by Heckman (1980) . The upper branch is identified 
with galaxies hosting Type 2 (Seyfert) AGN. 

Originally thought to be low-luminosity AGN, the na- 
ture of the LINER ionisation mechanism has since been 
debated, and it is becoming increasingly clear that the 
majority of galaxies with LINER spectra are unlikely to 

21 Both PPXF and GANDALF are available for download 
at http: //www. strw . leidenuniv .nl/saur on/ software .html. 

A version of GANDALF designed for use on SDSS spec- 
tra can be obtained at http://star-www.herts.ac.uk/ 
~sarzi/PaperV_nutshell/PaperV_nutshell .html 



be low-luminosity AGN. The physical footprint of the 
line emitting area in nearby LINERs has been deter- 
mined to be exten ded on kpc sc ales fe.g.. iKeell 119831 : 
[Phillips et al. 1986; Sarzi et al. 11200 6) and thus is unlikely 
be photoionized by a nuclear source. A detailed analysis 
of the distribution of this ext ended LINER era ission in 
nearby SAURON galaxies bv ISarzi etall ()2009D demon- 
strates that the majority of LINERs are incompatible 
with an ionization source in the nucleus, even in the 
case of galaxies known to host (X-ray or radio-detected) 
AGN, but instead are powered by some component^ of 
the underlying old stellar population. ISarzi et al.l p009( ) 
further show that this extended LINER emission is suf- 
ficiently strong to be detect ed in SDSS spectra of more 
distant z 0.05 galaxies. iCid Fernandes et"al] ([2009') 
have recently similarly argued that many of the SDSS 
objects in the LINER box are "retired" galaxies powered 
by old stellar populations. 

While it is certainly the case that some X-ray 
detected AGN exhibit nuc lear LINER spectra (e.g., 
[Gonzalez-Martin et all2009D . the presence of LINER- like 
line ratios over the large physical footprint of the SDSS 
spectroscopic fiber does not imply the presence of a low- 
luminosity AGN. LINERs are thus potentially highly di- 
verse and the majority of them are not low-luminosity 
AGN. Even if a low- luminosity AGN were present, the 
L[0 III] emission sampled by the spectroscopic fiber can- 
not disentangle the relative contributions. We therefore 
decided to eliminate them from our sample and to retain 
only the definite AGN for our analysis. 

Since galaxies exhibiting LINER spectra are different 
from those hosting Seyfert AGN, we verify that the inclu- 
sion of the small number of LINERs that are above the 
luminosity limit of L[Oiii] ^ lO""^ ergs~^established in 
the next Section does not significantly change our results, 
as most LINERs have substantially lower [O ill] A5007 
luminosities and Eddington ratios compared to Seyfert 
AGN and are thus not sites of substantial black hole 
growth (see Appendix [B] for details). 

The [O i]/Ha diagram allows for the best separation of 
AGN and LINERs, so we use it for those spectra where 
[Oi] A6300 is detected with S/N > 3. In cases where 
[Oi] A6300 is not detected, but [Sii] AA6717,6731 is, we 
use the [Sii]/Ha diagram (Figure [21 cen ter). In both 
cases, we use the AGN-LINER division of iKewlev eFall 
(|2006l) . If neither [Oi] A6300 nor [Six] AA6717,6731 is 
detected, but [Nil] A6584 is, we make use of the original 
[Nil] /Ha diag ram, where we us e the e mpirical division 
line defined bv ISchawinski et all ([2007D . 

Out of a parent sample of 47,675 galaxies, we thereby 
select a sample of 942 narrow line AGN (Table [SJ. This 
corresponds to a lower limit on the AGN fraction of 2%, 
as we do not include broad-line AGN, highly obscured 
AGN or LINERs. As we show later in this paper, the 
AGN fraction is highly dependent on position in the 
color-mass diagram and ranges from effectively zero to 
- 10%. 

3.3. Assessing the Completeness of Emission Line 
Selected AGN 

Since we wish to study the properties of AGN 
host galaxies as a class, we need to understand 
how effective our AGN selection technique is and 
to what degree it is biased. Emission line selec- 
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Fig. 2. — The emission line diagnostic diagrams used for the selection of AGN from narrow-line fluxes. We use three different diagrams, 
all of which use the ratio of [Olll]/H/3 {vertical axis). From left to right, the horizontal axes are the ratios of [NllJ/Ho, [Sn]/Hci and 
[Ol]/Ho. In each panel, we plot those objects that have S/N > 3 in each of the four lines required for the diagram. The gray shading 
represent the entire galaxy sample, while the green points are the AGN sample selected using these diagrams. In the [Nil] /Ha diagram 
(left), the dashed line shows t he empirical demar cation line between between purely star-forming galaxies and the composite region of 
the diagram determined by K auffmann et al.l ((2003a, Ka03). The solid line, labelled KeOl, is the theoretical extreme starburst line of 
[Kewley et al. (2001) beyond which the dominant source of ionization cannot be due to OB stars and must be due to a mechanism other 
than star formation (AGN, evolved stellar populations or shocks). The solid line labelled S07 is the empirical AGN-LINER distinction 
defined by Schawinski et al. (2007). The demarcation lines in the [Sll]/Ha and [Ol]/Ho diagrams are taken from Kewley et al. (2006). 
We select AGN using these diagrams in the following way: if [Ol] A6300 is detected with S/N > 3, then we use the [Oi]/Hq diagram. If 
[Ol] A6300 is not detected, but [Sll] AA6717, 6731 is, then we use the [Sii]/Hq diagram. If neither [Ol] A6300 nor [Sll] AA6717, 6731 are 
detected, but [Nil] A6584 is, we use the [Nil] /Ho diagram. 




Fig. 3. — [Nll]/Ho emission line diagnostic diagrams for all galaxies (shaded) and simulated composite AGN plus star forming galaxies 
(red points). The sim ulate d AGN luminosity increases from 10'^^ ergs~^in the left panel, to 10*" ergs^'^in the center and lO*'^ ergs~^on 
the right. See Section 13.31 for details. For low values of L[Olll] (left), the nebular emission from star formation can overwhelm the AGN 
signature, and such low-luminosity AGN may be missed. At an AGN L[Olll]= lO"" ergs"!, 12.8% of AGN move below the empirical 
starburst line (dashed), while 54.4% move below the extreme starburst line (solid). The result of this test shows that at high luminosities, 
AGN are efficiently selected using emission line diagrams. However, at lower luminosities, star formation may overwhelm the AGN narrow 
line emission. 



tion and characterisation of AGN has been success- 
fully e mployed in many wo r ks (e .g. . [Baldwin et ahl 
198lt IVeilleux fc Osterbroca [l98^ : IHo et al.l 119951 
1997a|lbllcl: IKauffmami et al.l I2003at [Sch awins ki et al.l 



2007t IConstantin et al.l 120081: iKewlev et al.ll200Q ). Deep 
Chandra studies of such emission-line selected AGN by 



iHo et all ()2001| ) show that the majority of such AGN also 
exhibit a nuclear X-ray source. For emission line selec- 
tion of AGN to be reliable and reasonably complete, two 
assumptions must be true: 

1. the AGN must always excite a narrow- line region, 
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-3-2-1012-3-2-10123 
Eddington parameter log,o Lpm/IVlBH 

Fig. 4. — Top left: The distribution of L[Oin] of our AGN sample. Top right: The distribution of stellar velocity dispersions of our AGN 
sample. Bottom: The distribution of the Eddington parameter L[Oin]/AfBHi combining the values from the plots on the top row via the 
A^BH — f relation. We indicate the approximate value of the Eddington limit in each panel. In each of the panels, we show the distribution 
of the entire AGN population in the top-left, while the remaining panels show the L[Olll] distribution split by host morphology. 



TABLE 2 

Typical [O in] A5007-luminosities, Eddington ratios and black hole masses of AGN by host galaxy morphology 



Host galaxy 


Number 


Percentage of 


mean 


median 


mean 


median 


mean 


median 


morphology 




host galaxies 


log L[Oiii] 


log L[0 III] 


log L[Oni]/o-'* 


log L[0 iii]/(t4 


log Mbh 


log Mbh 








( ergs-l) 


( ergs-l) 








(Mq) 


Late-types 


402 


42.68% 


40.40 


40.39 


0.295 


0.306 


6.52 


6.63 


Indeterminate-types 


431 


45.75% 


40.40 


40.39 


0.456 


0.442 


6.52 


6.63 


Early-types 


109 


11.57% 


40.27 


40.30 


0.326 


0.397 


6.36 


6.45 


All 


942 


100% 


40.38 


40.38 


0.372 


0.370 


6.43 


6.53 



and 

2. the emission lines from the narrow-line region must 
never be overwhelmed by lines excited by other pro- 
cesses, in particular star formation. 

IVIany recent studies of AGN host galaxies have noted 
that AGN tend to concentrate in the green valley on the 
color- mass diagram, i.e., at intermediate optical colors, 
and that there are very few AGN in the blue cloud and 



the red sequence. This raises the possibility that AGN 
in star-forming blue cloud galaxies could be overwhelmed 
by emission from star formation. The lack of relatively 
luminous AGN (Lx ^ lO'*-^"^^ ergs~^) in the blue cloud 
has be en recently been established by ISchawinski et al.) 
(120093) using hard X-ray-selected AGN, demonstrating 
that, at these luminosities, the absence of AGN is not a 
selection effect. 
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But could there still be a population of low-luminosity 
emission line AGN in the blue cloud whose emission 
lines are overwhelmed by star formation? For substan- 
tially lower lumi n osity A GN than those studied by by 
ISchawinski et all (|2009cl) . any X-ray detection in blue 
cloud star-forming galaxies is likely ambiguous in the 
0.2-10 keV window accessed by XMM-Newton and Chan- 
dra as this c ould be due to st ar formation and X-ray 
binaries (e.g.. iPersic et al.ll2004D . Can we do any better 
with emission line diagnostics? 

In order to test whether emission line AGN similar 
to those we select on BPT diagrams might be hidden 
in blue cloud star-forming galaxies, we perform a sim- 
ple empirical experiment. We pair a randomly selected 
star-forming galaxy from our galaxy sample with another 
randomly selected AGN. We rescale the AGN emission 
line fluxes to a specific i[Oiii] to simulate the effect of 
adding a AGN of a certain luminosity. For example, we 
compute the composite [N ii] A6584 line from the line of a 
random star-forming galaxy and the line from a random 
AGN as: 

[N ii]„,odci = [N ii]sF + f ) [N ii]agn (1) 

V -t^AGN / 

where Lg is the specific LfOiii] to which the AGN con- 
tribution to [N ii] A6584 is scaled. 

We perform this pairing of star-forming and AGN host 
galaxies 1,000 times scaling the AGN contribution to Lg 
= 10^9, 10-*° and lO"*! ergs"! in i[Oiii], to cover the 
range of star formation rate^ and line ratios of blue 
cloud star-forming galaxies and of AGP'O. We show the 
results of this exercise in Figure [3l where each panel 
shows a [Nil] /Ha diagram with the entire galaxy sam- 
ple as shades. On top of each, we plot the 1,000 com- 
posite objects scaled to the three specific L[Oiii] values. 
This addition of an AGN component moves the posi- 
tion of a star-forming object off the star-forming locus, 
into the composite region and finally onto the AGN re- 
gion. For added Sefyerts with L[Oiii] = 10*^^ ergs~^, we 
find th at only 1.4% remain below the iKaufFmann et al.l 
()2003aD empirical starburst line. However, for L[Olii] 
= 10"*° ergs~^and = 10^^ ergs~^, this fraction increases 
to 12.8% and 56.6% respectivel y. Similarly, th e fractions 
of AGN that move below the iKewlev et all (|2001[ ) ex- 
treme starburst line into the composite region are 85.5%, 
54.4% and 12.9% respectively. 

From this experiment, we conclude that towards the 
higher luminosities in our sample (~ 10^^ ergs^^), our 
sample is complete even in the blue cloud, and that there 
is no substantial population of such blue-cloud AGN. 
For lower luminosities, L[Oiii] < lO'*" ergs~^, this is 
no longer true, and there may potentially be AGN in 
the blue cloud whose narrow-line region is overwhelmed 
by current star formation. If there is such a popu- 
lation of AGN, they cannot be detected via emission 
line diagnostics, just as selection in the X-rays would 

The definition of "star forming" depends on the detection of 
emission fines with sufficient signal-to-noise and does not directly 
correspond to a star formation rate. 95% of star-forming galaxies 
have star formation rates greater than ~ 0.2MQyr~^ . 

Each pairing of an AGN with a star-forming galaxy is random, 
but we only have 942 AGN, so some of the AGN are used more 
than once. 



be challenging. However, a recent search for for broad- 
line A GN hidden in star-forming galaxies by iMao et all 
(|2009D found only three such objects in a parent sample 
of over 3,200 objects, establishing a lower limit on the 
very low-luminosity AGN fraction in the blue cloud of 
~ 0.1%, which suggests that the true abundance of such 
objects is very low, though broad lines may be similarly 
diluted by host galaxy light. 

3.4. Demographics of Black Hole Masses and Accretion 
Efficiencies 

We diagnose the accretion state of the black hole us- 
ing the [Qui] A5007 line. Observed correlations be- 
tween L[Oiii] and other indicators of AGN luminos- 
ity imply that i[Olll] is a reasonably good indicator 
of the c urrent accretion rate, though with considerable 
scatter (Heckman et al.' '2004) . More recent work by 
Diam ond-Stani c et al. (2009) comparing a variety of in- 
dicators of the intrinsic bolometric output of local AGN 
do show that L[Oiii] may not be an isotropic indicator 
and that Type 2 AGN - such as those in our sample - 
are intrinsically more obscured, and therefore more lu- 
minous, than pre viously thought, ev en b ased on hard X- 
ray o bservations (jRigbv et al.ll2009l : cf., ILamastra et al.l 
I2009D . Given these caveats, L[Oiii] will be used cau- 
tiously as a probe of black hole accretion in this work, 
and we restrict its use to statements about populations 
and trends rather than about individual objects. L[0 ill] 
is the only indicator of the current accretion rate in low- 
luminosity AGN available for large SDSS samples. 

We infer supermassive black hole masses indirectly 
from the stellar velocity dispersion at the effective 
radiu s via the Mrh — a rel a tion (iGebhardt et al] 
12000^, 'Ferrarese & MerrittI 120001 : iTremaine et al.l 120021 : 
H aring fc Rix 2004i') usin g the slope and normalisation 
of ITremaine et al.l (|2002l l . While the resolution of the 
SDSS spectra and the templates used is approximately 60 
kms~^, we are able to use GANDALF to probe slightly 
below this resolution limit. Whenever we quote black 
hole masses, we include all objects with velocity disper- 
sion measurements down to 40 kms~^, but caution that 
in the range from 40 - 50 kms~^(log Mbh= 5.3 — 5.7), 
the quoted black hole masses are increasingly uncertain. 
The median error in this lowest velocity dispersion range 
is ^ 10 kms"""^, corresponding to an error of ^ 0.4 dex 
in Mbh- 

To estimate the Eddington ratio of the AGN, we 
adopt the measure £[ ill] /Mrh as a rough proxy 
(jKewlev et al.l [20061) . iKauffmann fc HeckmanI (|2009f ) 
normalize the Eddington parameter to argue that the 
Eddington limit corresponds to log(L[0 iii]/Mbh ) = 
1.7 dex. Our AGN sample contains 33 (3.5%) objects 
with logL[0 iii]/A/bh > 1-7 dex, though the large scat- 
ter in L[0 iii]/A/bh m akes it likely that none o f our A GN 
are super-Eddington. IKauffmann fc HeckmanI (|2009f ) ar- 
gue that the bolometric correction for the extinction- 
corrected [O III] A5007 luminosity is a factor of 300 - 600, 
so the bolometric luminosities of the AGN in our sample 
are between 2.5 and 2.8 dex higher than the [O ill] A5007 
luminosities quoted throughout. 

We summarize some of the key AGN properties, in- 
cluding the AGN fraction in each morphology class, in 
Table |2l and plot the distribution of the stellar velocity 
dispersion, L[Oiii] and L[Oiii]/AfBH for both the entire 
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AGN sample, and for each morphology class, in Figure 
m We note that the majority of the AGN in our sample 
(~ 90%) reside in indeterminate- and late-type hosts. 

4. TWO DISTINCT MODES OF BLACK HOLE GROWTH IN 
EARLY- AND LATE-TYPE AGN HOST GALAXIES 

With the various observed and derived quantities as- 
sembled, we can now carry out a systematic investigation 
of the nature of AGN host galaxies and thus assess the 
role of AGN in the ongoing evolution of the low-redshift 
galaxy population. In Figure [SJ we show the location 
of AGN host galaxies on a color-mass diagram of m — r 
color vs. stellar mass as calculated in Section 12.41 In 
each panel, we show the entire normal galaxy popula- 
tion as contours. On top of the normal galaxy distribu- 
tion, we overplot the 942 AGN host galaxies as individ- 
ual green points to show their absolute distribution in 
the color-mass diagram. The top left panel contains the 
entire sample regardless of morphology. Going clockwise 
around the Figure, the top right panel shows only early- 
type galaxies, both for the normal galaxy population and 
the AGN host galaxies. Below in the bottom-right, we 
show only the late-types and finally in the bottom left, the 
indeterminate-types that did not receive a > 80% classi- 
fication. In the panels split by morphology, we overplot 
the contours of the entire population from the top left 
panel as gray dotted contours to aid the eye in identify- 
ing the underlying blue cloud and red sequence. 

The absolute distribution of AGN host galaxies on the 
color-mass diagram conveys the range of observed host 
properties, but can obscure the equally significant rela- 
tive incidence of AGN in various areas on the color-mass 
diagram. Previous studies have only used this absolute 
distribution to interpret the role of AGN in galaxy evo- 
lution - and then only for the whole galaxy population, 
not for the two main morphology classes separately. 

The AGN fraction on the color-mass diagram (Fig- 
ure [ni is more enlightening than the absolute distribu- 
tion because in any area on the color-mass diagram, the 
AGN fraction is a proxy for the duty cycle of AGN and 
therefore for the importance of AGN for the evolution 
of galaxies of specific masses, colors and morphologies. 
The duty cycle is the fraction of the time that a galaxy 
spends as a AGN in a specific population. The AGN 
fraction is a proxy of this because galaxies do not move 
rapidly on this diagram relative to the life time of AGN 
phases. Typical AGN lifetimes are estimated to be on the 
order of 10^ — 10® years, approaching 10^ years only for 
the most massive black holes which, in any case, a.re not 
active in the local Universe (iMartini fc Weinberg 120011 : 
iMarconi et"al] 120041: iMartinil |2004[ ). The movement of 
galaxies on the color-mass diagram on the other hand 
is slow: mass-doubling timescales du e to star formation 
are on the order of 10^ — 10^° years (jBrinchmann et ahl 
I200I INoeske et al.l 12001 while major m ergers are rare 
in the local Universe (jParg et al.ll2009b[ ). Stellar evolu- 
tion limits the pace of movement in the color direction. 
The most rapid color evolution possible is that of the 
instantaneous shutdown of star formation resulting in a 
movement from blue to red in ~ IGyr. Thus, the frac- 
tion of AGN in any area on the color-mass diagram is 
probing the fraction of time that AGN are switched on 
in that population, and therefore it is a measure of the 
duty-cycle. 



We present the AGN fraction for our sample in Fig- 
ure [6l which has a very similar layout as Figure [5) We 
do not fill the general galaxy population contours, but 
do leave the gray dotted contours indicating the whole 
galaxy population to guide the eye. We overplot the 
fraction of AGN host galaxies, with increasingly dark 
shades of green indicating a higher AGN fraction. It is 
immediately apparent from Figures [5] and |6] that the ab- 
solute and relative distribution of AGN host galaxies are 
substantially different, and that this difference is a very 
strong function of morphology. 

4.1. The Importance of AGN Host Galaxy Morphology 

Figures [5] and [6] suggest that the early-type and late- 
type host galaxies of AGN are fundamentally different 
populations, and that analysing the entire AGN host 
galaxy population as a whole obscures these differences. 
There is no 'typical' AGN host galaxy. 

We note at this point that, out of the entire sample 
of 942 AGN host galaxies, only 12 (1.3%) were classified 
as ongoing major mergers with serious disturbances by 
iDarg et "aD (|2009b( ): the vast majority of emission-line 
selected AGN in the local Universe are not hosted by 
galaxies undergoing a major merger. 

From their distribution in Figure [51 we see that the 
AGN host galaxies span nearly the entire observed 
range of the general galaxy population (in this sample) 
in stellar mass. This is not the case for color; many 
works have already noted that AGN host galaxies 
preferentially reside in the green valley at intermediate 
host galaxy colors and that they are le ss common in the 
blue cloud and the red sequence (e.g.. iSchawinski et al ' 
1 20071: [Nandra et al. 2007; Wild et all 120071: iSalim et al 
pot'; Consta ntin et all i2008: 'Martin et al.' '20071: 
| Georgakakis et ^l . | 20 081: Silverma n et al. 2002 
ISchawinski et al.ll2009cl : iTreister'et al.ir2q09bt ). 

The concentration of AGN host galaxies in the green 
valley is strongly dependent on whether we are consider- 
ing the relative or absolute AGN distribution. Figure [5] 
shows that a substantial number of AGN lie on the red se- 
quence. However, when we consider the AGN fraction in 
Figure [6l this population is put into context - compared 
to the total number of galaxies on the red sequence, the 
AGN population is lower by at least a factor of 10-20 
than in the peak regions, and therefore so is the duty 
cycle of these AGN. Hence, the importance of AGN to 
galaxy evolution on the red sequence is likely very small: 
we note that the population of AGN host galaxies at 
the massive end of the red sequence is almost entirely of 
indeterminate- and late-type morphology; massive early- 
type hosts are rare. 

The AGN fraction is strongly concentrated below the 
red sequence and at relatively high host stellar masses of 
M* ~ 10^"'^ Mq. Furthermore, the AGN fraction rapidly 
drops towards lower stellar masses and is very low to- 
wards the lower mass range. 

The absence of emission line AGN host galaxies on the 
red sequence is significant, as the emission line signature 
of even a very low-luminosity AGN superimposed on a 
passive stellar population would be detected. Interpret- 
ing the absence of AGN host galaxies from the blue cloud 
is more complex: As discussed in Section 13.31 their ab- 
sence is not significant for AGN at the lower end of the 
L[Oiii] of AGN in our sample (see Figure lU, because 
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Fig. 5. — The distribution of AGN host galaxies on the u — r color-mass diagram {green dots). Top left: galaxies of all morphologies; 
top right: early-type galaxies; bottom left: indeterminate-type galaxies; bottom right: late-type galaxies. The solid, shaded contours in 
each case show the galaxy population, on top of which we plot the individual AGN as green points. The large points are those AGN with 
L[Olll]> 10*^ ergs~^. In the panels for each specific morphology class, we only plot AGN host galaxies in that specific class. Furthermore, 
to guide the eye, we overplot the contours of all galaxies from the top left panel as gray dotted contours. The contour levels in all panels 
are linear and represent the same contour levels. From this Figure, we see that AGN host galaxies preferentially have green host galaxy 
colors and a range of stellar masses - but that this range is clearly morphology-dependent. That is, AGN host galaxies are a very particular 
sub-set of normal galaxies, and the way in which they are separate depends on morphology. 



the signature of the AGN narrow-hne region could be 
overwhelmed by star formation. At higher luminosities 
of > 10''^ ergs~^, the absence of blue cloud AGN is 
highly significant, as we would have been able to detect 
them. At even higher lu minosities of Lx ^ 10^'^ ergs~^, 
iSchawinski et al.l ( 2009cl) have established the complete 
absence of AGN in the blue cloud using hard X-ray ob- 
servations by the Swift satellite. 

What happens when we split the AGN population by 
morphology? In the remaining panels of Figure [S] and 
[6l a dramatically different picture emerges wherein the 
early- and late-type AGN host galaxies are shown to be 
significantly different populations. The absolute distri- 
bution of AGN host galaxies in Figure [5] already shows 
that the typical host stellar mass early-type AGN hosts 
is lower than that of late- type AGN hosts. 

Turning to the AGN fraction and therefore to where 
the AGN duty cycle is highest, this picture becomes even 
more extreme. The AGN fraction peaks at a stellar mass 



of ~ IO^^Mq in the general population, and ranges 
from 10^'^ — 10 Mq. Once we split by morphology 
in Figure [51 the difference between early- and late-type 
hosts becomes even more striking than for the absolute 
distribution: while for the late-types the AGN fraction 
still peaks strongly at 10^^ Mq, the AGN fraction in the 
early- type population peaks at M* ~ lO^^M©, an or- 
der of magnitude lower than in both the late-types and 
and the general population. This division has not been 
apparent in previous works due to a lack reliable morpho- 
logical classification for very large numbers of galaxies. 
Because the early- type AGN make up only 11% of all 
AGN (see Table [2|) , their properties are overwhelmed by 
the late-types in a general AGN host galaxies sample. 

The most striking aspect of the AGN fraction in Fig- 
ure [S] is that in the early- type population, the AGN duty 
cycle is highest in an extreme population of early-types 
that do not reside on the red sequence, but rather below 
it at the low-mass end and above the bulk of the blue 
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Fig. 6. — The distribution of the fraction of galaxies that host AGN on the color-mass diagram. The layout of this figure is similar to 
that in Figure[5] though for clarity, we do not shade the contours of the galaxy population. The filled contours represent the AGN fraction, 
which has been calculated in cells of size 0.175 dex in stellar mass and 0.125 mag in u — r color. We consider only cells that contain at least 
50 objects to minimize noise. The legend in the top left panel relates the shading level to the AGN fraction. As we argue in Section |4l the 
AGN fraction is a proxy for the AGN duty cycle. This Figure reveals in in stark terms that those galaxy populations that only very specific 
sub-populations of galaxies have a high AGN duty cycle: in the early-type population, it is only low-mass galaxies in the green valley that 
have a high AGN duty cycle. Early-type galaxies on the red sequence, regardless of mass, do not have an appreciable duty cycle. In the 
late-type population, it is massive galaxies that have a AGN high duty cycle. The duty cycle peaks strongly in a 'sweet spot' in the green 
valley. The high mass of late-type AGN host galaxies is somewhat misleading, as they have similar black hole masses, and therefore bulge 
masses, to the early-type AGN hosts. Their higher stellar mass is therefore due to a disk, rather than a bulge (see Table [2j. Example 
images are shown in Figure [TOl The Milky Way galaxy most likely resides in the 'sweet spot' for the AGN duty cycle in late-types. 



cloud, i.e., in the green valley. Similarly, the AGN duty 
cycle for late-type galaxies is highest in objects that are 
more massive and redder than the typical late-type pop- 
ulation. These red late-types will include a substantial 
fraction of inclined spirals with intrinsically blue colours, 
but dimmed and reddened by dust. If the detected AGN 
fraction is independent of host galaxy inclination, which 
is expected if the bulk of the nuclear obscuration comes 
from the central regions orientated randomly with re- 
spect to the galaxy, then these contaminants to the red 
part of the late-type colour-mass diagram from the part 
with much lower AGN fractions will make this conclusion 
even stronger. 

The indeterminate-type galaxies exhibit an AGN frac- 
tion distribution similar to the late-types, though their 
AGN fraction does not extend into the higher mass red 
sequence. This may be explained by the fact that the 



incidence of indeterminate-type galaxies in the red se- 
quence is much higher than that of late-types. We there- 
fore tentatively conclude that the indeterminate-type 
AGN, which are 45% of the AGN population after all, 
appear to behave more similar to the late-types than the 
early-types. 

When we make a cut in Eddington parameter of 
L[0 iii]/Mbh> 0.7, corresponding to an Eddington ratio 
of L/L^dd 0.1, an even more extreme picture occurs. 
By making a cut in Eddington ratio, we avoid any bias 
towards more massive black holes which radiate at higher 
luminosities at fixed L/L^dd than their low- mass coun- 
terparts. Figure [7] is identical to Figure [6] except for the 
cut in L/L^dd- This reveals that the duty cycle of high 
Eddington AGN is significant in only one population of 
low-redshift galaxies: low-mass early-type galaxies in the 
green valley. 
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Fig. 7. — The distribution of the fraetion of AGN host galaxies on the color-mass diagram for AGNthat have an Eddington parameter 
L[0 lll]/AfBH> 0.7 which corresponds to an Eddington ratio of > 0.1. The layout is identical to Figure[6l This Figure shows that the duty 
cycle oi^ high Eddington AGN is very low in all galaxy populations except for low mass early-types in the green valley. The prominent locus 
of high duty cycle in massive late-types in Figure |6] is therefore not dominated by high Eddington accretion. 



We conclude that the early- and late-type host 
galaxies of AGN are fundamentally different popula- 
tions: for late-type hosts, AGN preferentially reside in 
massive (IO^^Mq) galaxies either with green colors or 
on the red sequence. In contrast, early-type AGN tend 
to have moderate mass (10^'^ Mq ) and strongly cluster in 
the green valley in between the blue cloud and the red se- 
quence. We further showed that the duty cycle for high 
Eddington AGN is high only in green valley early-type 
galaxies. 

The AGN fraction on the color-mass diagram is a new 
and powerful tool to probe the duty cycle of AGN. Under 
the assumption that a high duty cycle AGN phase is more 
likely to affect the host galaxy evolution, the duty cycle 
is in turn a measure of the importance of AGN for galaxy 
formation for specific galaxy populations. 

4.2. Black Hole Demographics and Eddington Ratios in 
Early- and Late-type AGN Host Galaxies 

We have shown that the properties of early- and late- 
type AGN host galaxies are different. This in turn 
strongly suggests that the triggering and fuelling of the 



AGN, as well as their role in the host evolution, are differ- 
ent. We therefore continue our discussion with the black 
hole masses and accretion efhciencies of black holes in 
early- and late- type AGN host galaxies. 

Note that the characteristic black hole mass of the 
AGN is Mbh 10*^-5 Mq and that this value IS approxi- 
mately independent of the host galaxy morphology (Ta- 
ble [2]). The active black holes in late- types hosts are 
typically more massive than those in early-types by only 
^^0.2 dex, or less than a factor of 2. Combined with the 
fact that the typical accretion rate estimated by L[Oiii] 
is the same in early- and late-types, we find that the Ed- 
dington ratios as inferred via the Eddington parameter 
L[Oiii]/A/bh are similar, with the AGN in early-types 
having Eddington ratios that are only about a factor 1.5 
lower than those in late-types. Thus, the typical black 
hole mass growing today in early- and late-type galaxies 
is the same, with the same range of Eddington ratios. In 
addition, the host galaxy colors of the AGN host galaxies 
place them preferentially in the green valley, adding to 
an apparent picture where the conditions for black hole 
growth are approximately independent of host morphol- 
ogy- 
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Fig. 8. — The distribution of supermassive black hole masses Mbk for both normal galaxies and AGN host galaxies as a function of 
morphology. In each panel, we show histogram of A/bh for the normal galaxy population {gray, solid) and the AGN only (colored, dashed). 
Both here and in Figure [9] w e plot only objects where the measured velocity dispersion is greater than 40 kms~^ (corresponding to log{A/BH) 
~ 5.3: riYemaine et al.|[20M) . The mean and median black hole masses of both early- and late-type AGN host galaxies are very similar (see 
Table I2J, but when their distributions are compared to those of normal galaxies, large differences become apparent. Among the early- type 
galaxy population, it is preferentially those galaxies with lower black hole masses that are AGN, while the opposite is the case for the 
late- type galaxies. 



However, as we have shown in the previous section, 
the stellar masses of early- and late-type host galaxies 
are very different, so taking only population averages 
into consideration can obscure crucial variations between 
populations. Similarly, it is not just relevant which black 
holes are currently accreting, but also to what degree 
these growing black holes are different from those in their 
quiescent parent population, and crucially, how this dif- 
ference depends on host morphology. 

In Figure [51 we plot the distribution of inferred super- 
massive black hole masses for both the AGN (striped) 
and the normal galaxy population {soild), both for the 
entire sample, and split by host morphology. These his- 
tograms highlight that taking the mean, and ignoring 
the distribution of black hole masses in the parent pop- 
ulation masks intrinsically very different distributions of 
black hole growth. Even though the median black hole 
masses of active early- and late-type galaxies are compa- 
rable (Table [2]), a Kolmogorov-Smirnov test shows that 
they are drawn from different parent distributions at the 
99.7% significance level. 

The distribution of active black holes in late- type hosts 
skewed towards higher black hole mass, peaks at Mbh 
~ WMq and cuts off above Mbh ~ lO^M©. This con- 
trasts sharply with the distribution of black hole masses 
in all late- types galaxies, which decreases with increas- 
ing black hole mass. Plotting the fraction of active su- 
permassive black holes in late-type galaxies in Figure 
[9l we find a strongly increasing fraction from ~ 1% at 



lO^-^MQto ~ 8% at IQ'^-'^^Mq. At greater black hole 
masses, the fraction may begin to decline again, until we 
run out of number statistics. 

The black holes of early-type galaxies show a radically 
different picture (Figure |8]). The mass distribution of ac- 
tive black holes in early-type galaxies is skewed towards 
the lowest inferred masses in our sample. Plotting the 
active fraction as for the late-types in Figure [31 we find 
a very different behaviour. The fraction of active black 
holes in early-type galaxies is decreasing with increasing 
black hole mass, but we run out of galaxies to provide 
adequate number statistics. 

These two very different trends in early- and late-type 
hosts are unlikely to be selection effects. Since we are 
able to detect AGN in low-mass early-types down to host 
galaxies with stellar masses of 1Q'^'^Mq\ we do not see 
AGN in late-type hosts where the lowest-mass hosts are 
^ IO^-^Mq. In order to substantially change the trend 
of increasing AGN fraction with increasing black hole 
mass, we would have to have missed a very large number 
of AGN in low mass blue cloud late-type galaxies, more 
than ten times the number in our sample. To change 
the distribution to one skewed to low masses similar to 
the distribution of the early-types, an increase by a fac- 
tor of ~1000 would be required (Figure [H). At lumi- 
nosities above the ~ 10'"' ergs~^ completeness limit in 
Section [3.31 we are not missing any AGN in star-forming 
galaxies that might change this trend. Since the major- 
ity of the AGN contributing to the trends in Figure [SI 
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Fig. 9. — The fraction of supermassive black holes that are growing in the local Universe. We plot the black hole mass A/bh versus the 
fraction of galaxies that host an AGN. The errors on the fraction are Poisson errors. We plot the AGN fraction for the entire sample as 
well as for galaxies of specific morphology classes: (red) all; (green) early-types; (blue) late-types. This plot elaborates on the trend seen 
in Figure \8\ and shows the relative incidence of AGN in galaxies of early- and late-type morphology. The AGN fraction - and therefore 
the duty cycle - in early-type galaxies is heavily skewed toward low black hole masses. In late-type galaxies, the AGN fraction is skewed 
towards higher black hole masses, potentially with a drop at the highest masses. 



are above this limit (see Figure U), only lower-luminosity 
AGN might change these trends (see also|B]). Currently, 
there are no indications that such a large population of 
low-luminosity AGN exists (factor of ~1000 greater in 
number than observed) , but our results are an additional 
motivation to search for them. 

We restrict our AGN sample by Eddington ratio as 
we did in Section 14.11 by limiting the AGN sample to 
those AGN with L[Oiii]/Mbh> 0.7 (which corresponds 
to L/LEdd ~ 0.1) in order to test whether the trend we 
see in Figure [9] is due to a bias against low-mass black 
holes at a fixed Eddington ratio. Such a bias is unlikely 
to be the cause of the different trends we see for early- 
and late-type hosts as it would imply vastly different dis- 
tributions of Eddington ratios. Imposing this cut in Ed- 
dingto ratio to Figure |9] yields qualitatively the same 
trend, though the increase in the AGN fraction for the 
late- type hosts is less steep. 

In summary, dividing the local AGN host galaxy popu- 
lation by morphology reveals that that there are two very 
different modes of black hole growth: in early-types, it is 
preferentially the least massive black holes that are ac- 
tively accreting material in an A GN phase. In late-type 
galaxies, it is preferentially the most massive black holes 
that are active, with a potential decline at the very highest 
masses. Since the early- type AGN are only 11% of the 
population, this effect can easily be masked by studies of 
the AGN population as a whole. 

5. A NEW VIEW OF THE CO-EVOLUTION OF GALAXIES 
AND THEIR SUPERMASSIVE BLACK HOLES 



Despite the differences discussed in the previous sec- 
tion, the velocity dispersion - and therefore the bulge 
mass and black hole mass - of AGN host galaxies is in- 
dependent of morphology (Tableland Figure |4]). An- 
other common feature is the green host galaxy colors. 
This raises the possibility that these conditions are most 
favorable to fuelling a black hole, at least in low-redshift 
galaxies. The late- type AGN hosts are an order of mag- 
nitude more massive ^ IO^^Mq simply because of the 
presence of a massive disks. The bulges of late-types 
that host AGN appear to be very similar to the low- 
mass, ^ 10^° A/0, early-type galaxies that host AGN. 

However, as we have shown, the early- and late-type 
host galaxies of AGN are fundamentally different popu- 
lations both in terms of which galaxies host AGN, and 
also which black holes are growing. Only by comparing 
the AGN hosts to normal galaxies not just of the same 
mass and color, but also morphology, does it become 
clear that the two populations differ from each other in 
fundamental ways. 

What are the consequences of our findings for the role 
of AGN in galaxy evolution in the local Universe? In this 
section, we outline a new view of the impact of AGN on 
galaxy evolution and make the case that AGN perform 
very different roles in early- and late- type host galaxies. 

5.1. The Role of AGN in the Evolution of Early-type 

Galaxies 

As discussed in Section |4l while there are a number 
of AGN in early-type galaxies on the red sequence, the 
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Fig. 10. — Example SDSS gri images 80'''' x 80'''' in size of late- type AGN host galaxies with stellar masses above lO^M©. These images 
illustrate that massive late- type AGN hosts are not massive bulges with a small disk (SO or Sa) , but rather have small bulges and massive 
disks. 



fraction of early-type galaxies on the red sequence that 
host AGN is very small, i.e., the duty cycle of AGN in 
this population is small as well when compared to the 
peak duty cycle in the green valley. 

We see from Figure 15] that the duty cycle of AGN in 
early-type host galaxies is peaked in a well-defined band 
at intermediate green colors precisely in the green val- 
ley above the blue cloud and below the low mass end 
of the red sequence. If the early-type AGN hosts in the 
green valley are post-starburst objects, then there clearly 
is an ample supply of star-forming progenitor galaxies at 
the same mass in the blue cloud. These progenitors would 
evolve to the green valley if star formation is suppressed, 
and will continue to fade onto the low-mass end of the red 
sequence. Such pr ogenitors may already have early- type 
morphology (e.g.. iSchawinski et a l. 2009a,), or they may 
first transform their morpholo gy in a major merger, as 



simulations s e em to show fe.g.. [ Toomre fc Toomrg|1972l : 
iBarnesI 119921 : iBarnes fc HernquistI Such a ma- 

jor merger may also fuel the black hole, giving rise to 
the AGN phase observed at a later stage ( Sander s et al 



1 1988': 'Hop kins et 311120061) . Simulations bv iKavirai et al 
{2009) show that minor mergers with a gas-rich dwarf 
are not sufficient to move a red early-type progenitor to 
the optical blue cloud (though they may reach UV-blue 
colors) and so are not a good candidate for the trigger 
mechanism. 

Of course, it does not follow from the intermediate 
optical colors of green valley galaxies that they must 
be in transition from blue to red. There are multiple, 
very different star formation histories that can result in 
green colors. A slightly enhanced dust screen covering 
an actively star forming galaxy would suffice. More com- 
plex scenarios are also imaginable. That early-type AGN 
hosts gal axies genuinely are a tra nsition population was 
shown bv lSchawinski et al.l(j2007D : however, this does not 
imply that the AGN phase in the green valley is necessar- 
ily the cause for the shutdown of star formation and the 
transition from blue to red. Rather, the absence of active 
bl ack hole accretion in blu e cloud objects, established 
bv lSchawinski et al.l (|2009c|) . implies that there must be 
a significant time delay on the order of at least 100 
Myr, more likely more like 500 Myr, between whichever 



process suppresses star formation and the onset of sub- 
stantial AGN radiation. Furthermore, ISchawinski et al.) 
(|2009bf ) showed that the extremely rapid destruction of 
the molecular gas reservoir that fuels star formation co- 
incides with a very low-level of accretion preceding the 
AGN phase in so-called composite objectto and that this 
very low-luminosity phase ~ perhaps acc ompanied by a 
radiatively inefficient kinetic outfiow (e.g.. lNaravan fc Yil 
[1991 - may be responsible for the shutdown of star for- 
mation by removing the molecular gas reservoir. 

Regardless of whether any AGN phase is responsible 
for the shutdown of star formation, the duty cycle of 
AGN in the present day early-type galaxy population is 
highest in low-mass green valley objects in which star 
formation has already been shut down and which are in 
the process of moving to the red sequence. The duty 
cycle of AGN in massive, red sequence early-types on 
the other hand is very low. The majority of black hole 
growth in early- type galaxies at low redshift thus appears 
to be a post-starburst phenomenon, rather than part of 
a coeval starburst. Figure [7] furthermore illustrates that 
it is only this population of low redshift galaxies that has 
a substantial duty cycle of high Eddington accretion. 

5.2. The Role oj AGN in the Evolution oj Late-type 
Galaxies 

The late-type AGN host galaxies are genuine late-type 
spirals, not SO/Sa galaxies that consist of a large bulge 
and a small disk, even at the highest masses, where the 
duty cycle is highest. We show example images of mas- 
sive late-type AGN hosts in Figure [TUl While the host 
galaxies are quite massive, their bulges are not - this is 
apparent from their bulge velocity dispersions, and there- 
fore their modest black hole masses, which are compara- 
ble to those of active early-type galaxies; they therefore 
must have substantial disks which make up the remain- 
ing stellar mass (Figure [5] and Table [5]) . To put it an- 
other way, the typical black hole masses, bulge masses 
and Eddington ratios of early- and late-type AGN hosts 

^* Composite objects lie on the BPT diagram in betwreen the 
locus of star forming galaxies and the non-stellar (AGN & LINER) 
locus and are usually interpreted as a combination of star formation 
and very low-level AGN activity. 
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are similar, but the late- type hosts are substantially more 
massive due to the presence of a massive disk. 

At the same time, the distribution of Eddington pa- 
rameters (Figure 2] and Table [2|) does not provide a com- 
prehensive picture: comparing Figure [6] shows that the 
duty cycle of AGN - as we are able to detect them via 
emission line diagnostics - is high in high-mass, green 
late- types. However, Figure [7] shows that the duty cy- 
cle of high Eddington-rate accretion is low in those late- 
types, especially when compared to the green valley 
early-types. 

Given the distinct location of late-type AGN host 
galaxies on the color-mass diagram, what role do the 
AGN play in the evolution of late-type galaxies? We ar- 
gued previously that the early-type AGN represent the 
true green valley transition population. The majority of 
AGN hosted by late-type galaxies do not reside between 
the blue cloud and the low-mass end of the green valley 
and therefore cannot be part of the same blue-green-red 
transition. A minority of the late-type AGN hosts do 
however have sufficiently low masses to plausibly origi- 
nate from the blue cloud. The absence of a substantial 
duty cycle for high Eddington accretion in these objects 
further supports this, as the AGN of the early-type hosts 
do have such a high duty cycle. There is also an absence 
of AGN at the very highest masses > 10 ^^'^M© which 
may indicate a feedback process that has shut down ac- 
cretion in those objects, or this may be due to insufficient 
number statistics. 

We consider a number of scenarios that may account 
for the presence of AGN in late- types; any such sce- 
nario must account for three observational character- 
istics: (i) characteristic host mass of ~ IO^^Mq, (ii) 
late-type host galaxy morphology, and (iii) intermedi- 
ate/green host galaxy color. 

5.2.1. Major Mergers of two Disk Galaxies 

The first possible scenario for the origin of the late- 
type AGN hosts is the major merger of two disk 
galaxi es that results in the fuelhng of the blac k hole 
(e^'Springcl & Hcrnqiiil [2005t iHopkins et all 120051: 
[Robert son et al. 2006). We know from the visual in- 
spection of merging galaxies in the SDSS U niverse that 
disk- d isk mergers at ~ IO^^Mq do occur (|Darg et al.l 
l2009al |bl). However, there is a lack of blue progenitor 
galaxies at ^ IO^^Mq that might evolve toward the red 
sequence as star formation is suppressed (see Figure |6]). 
Both factors thus make it unlikely that late-type AGN 
hosts are the product of such a merger and so we con- 
clude that this scenario does not account for criterion 

Also, a major merger would likely have destroyed the 
disk that is clearly visible in the SDSS images (Figure 
nop . Simulations show that mergers involving particu- 
larly gas-r ich progenitors can result in a disk-dominated 
remnant (ISpringel fc Hernguistl l2005t [Robertson et al.l 
[20061: IHopkins et al.ll2009ir , but only for particular orien- 
tation of the angular momentum vectors. Major mergers 
are already rare, and this requirement further reduces 
that probability by a large factor. Moreover, such ex- 
tremely gas-rich, massive progenitor galaxies do not ex- 
ist in the low-redshift Universe (e.g. lKannappanll2004D 
and so major mergers are an unlikely channel for the 
triggering of AGN in low redshift late- type galaxies. 



5.2.2. Mixed-mergers of a Disk and a Sptieroid Galaxy 

One possible scenario for the origin of late-type AGN 
is a mixed merger between a disk- and a massive 
spheroid-dominated system, which is predicted to be 
very common by some models of galaxy formation (e.g., 
iKhochfar fc Burkertl [200l iKhochfar fc Silkl l2006f) . The 
mixed merger scenario alleviates one of the problems 
with major mergers, as there is an ample supply of mas- 
sive early-type, red progenitors available on the red se- 
quence. It also accounts to some extent for the green 
host galaxy colors (iii); while the red sequence progeni- 
tor does not have a ready supply of cold gas available for 
star formation, the late-type progenitor does. 

However, this scenario also faces some problems: why 
does the mixed merger result in green, rather than blue 
optical colors? One explanation might be that the late- 
type progenitor originates from the blue cloud proper 
and is substantially less massive (~ bringing in 

only a small amount of cold gas that can fuel only a minor 
starburst. Simulations of minor mergers between pas- 
sive early-t ype galaxies an d an en semble of small gas-rich 
galaxies bv iKavirai et al.l ()2009D show that such events 
can produce the intermediate u — r colors exhibited by 
the late-type AGN hosts. 

While a major mixed merger likely also destroys the 
disk in the late- type progenitor, a minor merger might 
not, thus satisfying criterion (ii). However, another prob- 
lem is that the bulges of the late-type AGN hosts are 
comparable in velocity dispersion to the significantly less 
massive early-type AGN hosts and therefore are predom- 
inantly disk-dominated systems; a mixed merger of this 
type seems unlikely to build a massive disk (similar to 
that of the Milky Way; see Section 15. 5p around a small 
bulge. Thus a contradiction arises: in order to account 
for the green colors, the early-type must be more massive 
than the late-type, but in order to ensure the survival of 
the disk, the late-type must be more massive than the 
early-type. 

Another possibility is that the AGN is regulating the 
strength of the starburst and preventing it from moving 
the host galaxy all the way to blue cloud colors, which is 
commensurate with the idea that AGN are more common 
in late- type galaxies with high black hole masses. 

5.2.3. Concluding Remarks concerning Late-type AGN Host 

Galaxies 

All the scenarios outlined here are problematic and 
none of them account for all the observed properties of 
late- type AGN host galaxies. Perhaps an AGN phase in 
late-type galaxies does not significantly affect the evo- 
lution of the host galaxy, and all we see is secular evo- 
lution. A slight increase in the amount of gas in the 
bulge increases the odds for stochastic accretion of gas , 
perhaps drive n by a bar (jHasan et al.lll993l: lWvsell2Q0l: 
iConibesI 120061) or more gene r ally b y disk-driven evolu- 
tion iKormendv fc KennicuttI (|2004D . Othe r possibilities 
include mass loss from evolved stars (e.g., IDavies et al.l 
|2007 : Hq 2008) , the infall of a minor dwarf galaxy or 
the tidal disrupti on of a star by the black hole (e.g., 
IGezari et al.l[20Ci6l) 

As we have discussed, it is unlikely that the late-type 
AGN hosts are involved in the same transition from the 
blue cloud to the red sequence that the early-types ap- 
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pear to be involved in. The role of AGN in the evolution 
of late-type galaxies must therefore be a different one 
from the early-types. A major difference between the 
early- and late-type AGN hosts that highlights the dif- 
ferent nature of black hole accretion is that the late- types 
do not have any sub-population that has a high duty cy- 
cle of high Eddington accretion (Figure [T]) . Since there 
are no major ongoing starbursts in late-type AGN hosts 
and the duty cycle for high Eddington accretion is low, 
perhaps all back hole growth has already taken place at 
high redshift and all we are seeing is a minor re-activation 
due to one of the processes mentioned above. 

Given that none of the scenario we outlined are satis- 
factory explanations, the co-evolution of late-type galax- 
ies and their black holes ~ unlike in the case of early-types 
- remains an outstanding issue. 

5.3. Parallels between Black Hole Growth in the 
Present-day and Early Universe 

The only part of the early-type galaxy population that 
is experiencing a high duty cycle of black hole growth at 
high efficiencies are low-mass post-starburst early-type 
galaxies transitioning from the blue cloud to the low- 
mass end of the red sequence. This population is con- 
tinuing the build-up of the red sequence, so it is natural 
to speculate that the process is a downsized version of 
the process that led to the formation and quenching of 
more massive early-type galaxies at high redshift which 
are known to be older and to have formed on shorter 
time scales. 

The late-type AGN are not associated with the same 
transition from blue cloud to the red sequence and their 
role in the evolution of their hosts is less clear. The ma- 
jority of the black hole growth in the local Universe is 
associated with late-type hosts (at least ^ 40%, and al- 
most 90%, if we include all indeterminate- types), and 
only ~11% is associated with early-types. Thus, the 
mode of black hole growth at work in late-type galax- 
ies is significantly more common in the local Universe 
than the one at work in early-types, and so most black 
hole growth is not associated with the blue-red transition 
seen in early-types. 

Recent observations o f X-ray-selected AG N host galax- 
ies in deep fields by iGeorgakakis et all (|2009f ) sug- 
gest that the fraction of AGN host galaxies with late- 
type morphology is only ~ 30% at z ~ 1, indicat- 
ing a decrease of the late-type fraction with increas- 
ing redshift, thou g h the AGN population studied by 
IGeorgakakis et aT] (|2009f ) is somewhat more luminous 
that in our sample. This decrease may suggest that the 
appearance of the mode of black hole growth at work in 
late-type galaxies is relatively recent, perhaps becoming 
apparent only around z ^ 1 when the assembly of the 
Hubble sequence begins to finalize. Another possibility 
is that this mode has always been active at a moderate 
level, but that the AGN host galaxy population at high 
redshift is dominated by the downsizing mode at work 
in the early-type population. The two modes at work 
in early- and late-type galaxies, respectively, change in 
importance with cosmic time, the late-type mode may 
account for the the observed peak at z ~ 0.7 in the num- 
ber density of moderate-luminosity X-ray selected AGN 
at z 0.7, while the mode at work in late-types is re- 
sponsible for the high redshift peak of luminous AGN 



(iGiacconi et all [200 It [Alexander et~alll2001l: lUeda et al.l 
12003 ; iGilh et al.ll2005D . As star formation and black hole 
growth moves to less massive early- type galaxies (e.g., 
iThomas et"aI1[2 005. 2009), the black hole growth mode in 
early-type galaxies also becomes less luminous and thus 
by < 1 becomes overwhelmed by the late-type mode. 

Perhaps the local Universe does not tell us much 
about the AGN-host galaxy connection because black 
hole growth rates and and star formation rates have 
declined so substantially since z ~ 1 except in a few 
rare systems more reminiscent of activity at high redshift 
(e.g. NGC 6240)? Further work on the morphology of 
AGN host galaxies at high redshift, in particular detailed 
restframe optical morphological classification, will be re- 
quired to investigate these suggestion. Near-IR imaging 
data from the new Wide Field Camera 3 on board Hubble 
could be used to reproduce this work at z ^ 2. 

5.4. What about Indeterminate-type AGN Host 
Galaxies ? 

Throughout our entire discussion, we have mostly con- 
sidered the early- and late-type samples as the two mor- 
phological extremes. However, both in terms of the entire 
galaxy sample and the AGN sample, the indeterminate- 
types are the most numerous class (45.3% of galaxies and 
44.5% of AGN; see Tables [T] and H]). Future imaging data, 
such as from PanSTARRS, VST0 or LSSTEH, classi- 
fied via a future Galaxy Zoo project would enable us to 
better separate them into early- and late-types. 

5.5. The Black Hole at the Genter of the Milky Way 
should be Active! 

Figure [6] shows that the duty cycle for AGN in late- 
type galaxies is very strongly peaked at a particular stel- 
lar mass and color. The Milky Way most likely occu- 
pies this sweet spot for late-type AGN. Estimates of 
the stellar mass of the Milky Way vary, but are around 
^ IO^^Mq, while the current star formation rate is 
approxim ately 3 Mpyr~^ (e.g.. iScalol 119861; iDiehl et "aD 
2006; Calz etti et al.l 120091 ). These parameters place the 
Milky Way right where the AGN fraction for late-type 
galaxies is highest, between 5 and 10% and therefore the 
duty cycle is the highest. However, this also implies that 
the black holes in galaxies like the Milky Way are not 
active ~ 92 — 95% of the time at the luminosity of the 
typical Seyfert in the SDSS Universe (~ 10^° ergs""'^, see 
Section [S3] and 

The black hole at the center of the Milky Way at Sagit- 
tarius A* (Sgr A*) is, apart from occasional weak flares, 
remarkably quiet with an extremely sub-E ddington ac- 
cretion rate (e.g., [MaiEQneIeFn]|200l, [2008|). However, 
there is evidence that Sgr A* was substantially more ac- 
tive in the very recent past. Hard X-ray observations 
with ASGA showed Fe Ka lines around Galactic Cen- 
ter emitted by the molecular cloud complex surrounding 
the black hole. Follow-up observations with more recent 

The Panoramic Survey Telescope & Rapid Response System; 
see http: //pan-Starrs . if a.hawaii . edu/public/. 

26 The V LT Survey Telescope ; 

http : / /www. eso . org/ sci/observing/policies/PublicSurveys/ 
sciencePublicSurveys .html#VST. 

2^ The Large Synoptic Survey Telescope; 
http: //www. Isst . org/lsst. 
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TABLE 3 

A CATALOGUE OF AGN HOST GALAXIES 



ID 


SDSS Object 


RA 


Dec 


GZ 


Redshift 


A^stollar 


cr 


u — r 


logL[Oni] 


logL[Oni]/cr4 




ID 


{J2000) 


{J2000) 


class^ 




Mq 




colour 


ergs^^ 




1 


587722983902544102 


14 17 04.4 


+00 30 28.8 


i 


0.025730 


10.01 


69.9 ± 4.91 


2.51 


39.53 


-0.181 


2 


587724197207277638 


01 31 01.1 


+ 13 03 15.7 


e 


0.020734 


10.80 


115.5 ± 3.23 


2.60 


39.33 


-1.308 


3 


587722984428929165 


12 40 53.6 


+01 00 30.1 


e 


0.022932 


10.76 


73.7 ± 2.92 


2.64 


39.86 


0.033 


4 


587724232647704698 


02 14 04.4 


+ 13 11 56.4 


i 


0.041659 


10.40 


106.4 ± 6.18 


2.16 


41.25 


0.853 


5 


587724233712205954 


00 47 30.3 


+ 15 41 49.5 


1 


0.031463 


11.21 


128.8 ± 3.08 


2.59 


39.94 


-0.819 


6 


587724231569309775 


01 30 06.2 


+ 13 17 02.1 


1 


0.037992 


10.16 


36.6 ± 10.60 


1.90 


40.69 


2.153 


7 


587724234257924238 


02 11 21.8 


+ 14 30 15.5 


1 


0.041690 


10.56 


81.4 ± 5.73 


2.44 


40.72 


0.876 



Note. — Table[3]is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding 
its form and content. 

^ Galaxy Zoo morphology: c - carly-typc; 1 - late-type, i - indeterminate-type. See Section l2.3l 



hard X-ray facilities indicate that Sgr A* was orders of 
magnitude, perhaps 10^ times more luminous in the last 
few hun dred years than i t is cu rrently. Suzaku observa- 
tions by iNobukawa et al.l (|2008[ ) detect a Ka luminosity 
consistent with Sgr A* having been at IQ^s-SQ ergs~^ a 
mere 300 yrs ago, while iRevnivtsev et al.l (|2004D used IN- 
TEGRAL data to argue that Sgr A* was at a 2-200 keV 
luminosity of L ~ 1.5 x 10^^ ergs^^ just 300 to 400 years 
ago. This is not quite the luminosity necessary to excite 
a narrow-line region sufficiently luminous to be detected 
by an SDSS fibre at z ^ 0.02, but these observations 
clearly imply that Sgr A* is active at slightly lower lumi- 
nosity on short timescales. Together with the duty cycle 
for Milky Way-like galaxies derived here, these observa- 
tions may give a first hint of the lifetime of AGN in such 
galaxies. A caveat to this statement is the fact that the 
luminosity of this very recent fiare is still below the AGN 
luminosity limit of our sample. 

Under what conditions would the black hole at the 
center of the Milky Way reach a sufficient luminosity 
so that it would be included in our AGN sample? To 
reach a typical luminosity of a few percent Eddington 
(~ few X lO''^ ergs~^ bolometric luminosity) requires 
a mass-accretion rate of only lQ~^MQyr~^ . This only 
a factor 100 higher than Bondi accretion rate inferred 
from simulations, and there are IQ'^Mc^, of molecular gas 
within 2 pc of the central black hole (jHerrnstein fc Hoi 
|2002|) . with large clouds at further distances. It would 
not necessarily be surprising if dynamical friction or a mi- 
nor merger, e.g., the swallowing of the Sagittarius dwarf 
galaxy drove some of these clouds in, temporarily re- 
igniting the black hole at the center of the Milky Way. 
All the ingredients thus appear to be on hand for turning 
our galaxy into a Seyfert AGN. Indeed, the mystery may 
rather be why is our central black hole so dark and why 
is gas accretion is so infrequent? 

The implication of this work is that the Milky Way is 
in a location on the color-mass diagram where late-type 
galaxies have the highest duty cycle for AGN phases. The 
mass of the Mil ky Way black hol e has been inferred to 
be 4.1 X 1O^M0 (jGhez et al.ll2008[ ) which places it in the 
typical range of black hole masses in late-type galaxies 
that are currently growing in the local Universe (Figure 
[5]and[Hl). Thus, 5-10% duty cycle for galaxies like the 
Milky Way places into context the current quiescence 
of Sgr A*, together with the evidence for accretion in 
the recent past and apparently favorable conditions for 



a resumption of accretion. 

6. SUMMARY 

In order to understand the co-evolution of galaxies and 
the supermassive black holes at their centers, we have 
investigated which galaxies are more likely to host ac- 
tive black holes and how they differ systematically from 
their normal counterparts, if at all. We use data from 
the Sloan Digital Sky Survey and visual classifications 
of morphology from the Galaxy Zoo project to analyze 
black hole growth in the nearby Universe {z < 0.05) and 
dissect the AGN host galaxy population by color, stellar 
mass and morphology. We show for the first time the im- 
portance of host galaxy morphology for black hole growth 
and relate it to our understanding of the way in which 
galaxies and their supermassive black holes co-evolve. In 
summary: 

1. We confirm that the selection of AGN via emission 
line diagnostic diagrams may miss low-luminosity 
AGN (L[Oiii] < 10^^° ergs"i) in star-forming 
galaxies. 

2. AGN host galaxies as a population have high stellar 
masses around IO^^Mq reside in the green valley, 
and have median black hole masses around Mbh 
- IO^-^Mq. 

3. When we divide the AGN host galaxy popula- 
tion into early- and late-types, key differences ap- 
pear. While both early- and late-type AGN host 
galaxies have similar typical black hole masses and 
Eddington ratios, their stellar masses are differ- 
ent, with the early- type hosts being significantly 
less massive (~ lO^^M©) than the late-type hosts 
(~ IO^^Mq). The late-type hosts nevertheless have 
the same characteristic bulge masses, and therefore 
black hole masses, as the early-types. The differ- 
ence in stellar mass is due to the presence of a mas- 
sive stellar disk in the late- types. 

4. We consider the fraction of galaxies hosting an 
AGN as a function of black hole mass. Dividing by 
morphology into early- and late- type galaxies, we 
find that in early-type galaxies, it is preferentially 
the galaxies with the least massive black holes that 
are active. In late-type galaxies, it is preferentially 
the most massive black holes that are active, with 
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a potential drop in the active fraction above (A/bh 
- IO^-^Mq). 

5. We estimate the duty cycle of AGN on the color- 
mass diagram, split by morphology. While some 
early-type galaxies lie on the red sequence that host 
AGN, their duty cycle is negligible. The duty cycle 
of AGN in early-type galaxies is strongly peaked 
in the green valley below the low-mass end of the 
red sequence and above the blue cloud. The duty 
cycle of AGN in late-types, on the other hand, is 
high in massive green and red late-types. When 
we impose a minimum Eddington ratio cut, we find 
that only the green valley early-type AGN still have 
a substantial duty cycle. No part of the late-type 
population has a high duty cycle of high-Eddington 
ratio accretion. 

6. We discuss the implications of these results for our 
understanding of the role of AGN in the evolution 
of galaxies. We conclude, in particular, that there 
are two fundamentally different modes of black hole 
growth at work in early- and late-type galaxies. We 
now have a good understanding of the role of AGN 
in early-type galaxies: The low-mass early-type 
AGN hosts are post-starburst objects moving to- 
wards the low-mass end of the red sequence. They 
may represent a low-mass, low-luminosity or down- 
sized version of the mode of evolution that was in- 
volved in the formation of more massive early-type 
galaxies at high redshift. The role of AGN in the 
evolution of the late-type galaxies is less clear and 
we offer some possible scenarios. The high host 
stellar masses and grand-design stellar disks makes 
it implausible that they are the product of recent 
major mergers moving from the blue cloud to the 
red sequence. 

7. Finally, the Milky Way likely resides in the 'sweet 
spot' on the color- mass diagram where the duty 
cycle of AGN is highest amongst late-type galaxies 
and that the AGN duty cycle for galaxies like the 
Milky Way is - 5 - 8%. 

We provide the entire AGN host galaxy sample used 
in this paper, including derived properties such as stellar 
masses, black hole masses, [O ill] luminosities and Ed- 
dington parameters, in Table [3l 
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Fig. 11. — The effect of AGN contribution to stellar mass estimates from optical photometric data. We show the results of a test of 
the effect of adding a quasar template (from Vandcn Berk ct al. 2003) to an ensemble of galaxy templates varying the relative V-band 
luminosity of the quasar and the galaxy template. This test shows that we are able to recover the stellar mass within 0.2 dex for galaxies at 
least ten times more luminous in V than the quasar. Stellar mass measurements of host galaxies of brighter quasars degrade as the quasar 
luminosity approaches the host galaxy. 



APPENDIX 

RELIABILITY OF STELLAR MASS MEASUREMENTS OF AGN HOST GALAXIES 

The reliability of galaxy parameters determined from s pectral energy distribution fitting techniques is limited by 
uncertainties in the input models and parameters (e.g., Maraston et al. 2006; Conroy ct al. 2009a, b). As we perform 
this stellar mass measurement not only for normal galaxies, but also for the host galaxies of obscured AGN, we 
need to assess the effect of he presence of a small amount of AGN continuum on stellar mass measurements. We 
generate a set of rnodel S EDs by combining a galaxy template, /a, galaxy with a quasar template /a, quasar 7 (taken from 
iVanden Berk et aLll2001l ) and normalising it by a constant a to achieve a specific ratio of iy_quasar/-^y,gaiaxy in order 
to explore a range of relative luminosities, specified by the constant, a: 



^A, model — /a, galaxy ~t~ Q^J'a, quasar 

(Al) 

The galaxy templates used are SSPs (single stellar populations) with ages of 100 Myr (young), 1 Gyr (intermediate) 
and 10 Gyr (old) and a 1 Gyr old exponentially declining CSP (composite stellar population). The results of this 
experiment are shown in Figure 1111 We find that for the older galaxy templates, the stellar mass estimates remain 
relatively unaffected up to a luminosity ratio of Ly,quasar/iy,gaiaxy 0.1. The stellar mass estimates for the younger 
templates begin to become uncertain at the level of ~ 0.2 dex when iy.quasar/iy,gaiaxy reaches ~ 0.01, which is likely 
the maximum amount of contamination that would be expected in our sample. 

Adding the quasar template to the galaxy is comparable to adding a very young stellar population. This explains 
the perhaps counterintuitive outcome of old stellar populations being relatively unaffected - the two-component code 
identifies the quasar as a small, very young stellar population with a very low mass-to-light ratio and isolates it from 
the bulk old stellar population. The SED fitting method has more trouble with galaxy templates where young stellar 
populations are genuinely present. In this case, the quasar template becomes entangled with the young population 
and so leads the code to identify a young stellar population of the wrong age and mass-fraction, resulting in a greater 
uncertainty in the measured stellar mass. This means that this effect is not a weakness of two-component codes, but 
a generic problem for all methods. 

This experiment implies that the uncertainty introduced by AGN is comparable to the other uncertainties that affect 
stellar mass measurements. The AGN in the relatively local volume probed by our sample are low-luminosity and 
obscured, so their expected AGN continuum contribution is small. As a consequence, the stellar mass estimates of 
AGN host galaxies are approximately as reliable as those of normal galaxies. 

ASSESSMENT OF THE IMPACT OF LOW-LUMINOSITY AGN AND LINERS 

Since LINER spectra are preferentially detected in massive red sequence early-type galaxies (e.g., iKewlev et al.l 
I2OO6I : iSchawin ski et al.l 12007*) . the exclusion of LINERs may bias our AGN host galaxies sample. If LINERs are 
AGN, t hen th eir im plied Eddington ratios inferred from _L[Oiii]/AfBH are systematically lower than those of Seyferts 
(|Kewlev et all [20061 ) . This implies that their [Oiii] A5007 luminosities are substantially lower for a given black hole 
mass, which in turn means that the tpyical LINER [O ill] A5007 luminosity is significantly lower than that of Seyferts. 
In fact, only 19% of all LINERs detected using our BPT diagrams have [Qui] A5007 luminosities greater than the 
completeness limit of ^ 10"'" ergs~^that we determine in Section l373l In order to test whether our removal of LINERs 
from our sample significantly affects our results, we perform the following experiment: we include all LINERs in the 
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Fig. 12. — Same are Figure|6l but for all objects classified as Seyfert or LINER with L[Oin]> lO'*" ergs~^. This Figure illustrates that 
the inclusion of the small fraction of LINERs that have high [O III] A5007 luminosities does not substantially change the results presented 
in this work. 

AGN sample, but remove all objects (whether classified Seyfert or LINER) with L[Oiii]< 10""^ ergs~^ (this removes 
^20% of the Seyfert sample as well). When we then redo the analysis of this paper using this new sample, we find 
that all the qualitative trends we describe remain unchanged, and that the quantitative differences are small. 



